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Abstract. We study the conductance through two types of graphene nanostructures: 
nanoribbon junctions in which the width changes from wide to narrow, and curved 
nanoribbons. In the wide-narrow structures, substantial reflection occurs from the wide-narrow 
interface, in contrast to the behavior of the much studied electron gas waveguides. In the 
curved nanoribbons, the conductance is very sensitive to details such as whether regions of 
a semiconducting armchair nanoribbon are included in the curved structure - such regions 
strongly suppress the conductance. Surprisingly, this suppression is not due to the band gap 
of the semiconducting nanoribbon, but is linked to the valley degree of freedom. Though we 
study these effects in the simplest contexts, they can be expected to occur for more complicated 
structures, and we show results for rings as well. We conclude that experience from electron 
gas waveguides does 7iot carry over to graphene nanostructures. The interior interfaces causing 
extra scattering result from the extra effective degrees of freedom of the graphene structure, 
namely the valley and sublattice pseudospins. 

PACS numbers: 73.63.Nm, 73.21.Hb, 73.23.Ad, 73.61.Wp 

1. Introduction 

There has been tremendous interest recently in investigating carbon-based nanoelectronics, 
first with carbon nanotubes HI HO and more recently with graphene p]. In that context, 
researchers have intensively studied graphene "nanoribbons" - infinite, straight strips of 
graphene of constant width - both theoretically |l5l|S|71[8l|9l[l0l[IT][I2[l3l[Il[I21[l6l 
[niiniEllEQiETllaalESIand experimentally ||24l|25]|2i|221|2a|29l[30l[3lll. Most of the 
theoretical effort has been focused on nanoribbons of essentially constant width. However, 
more functionality, beyond that of a mere wire, might be gained if one considers more general 
and realistic nanoribbons in which the width of the ribbon changes, it curves, or particular 
junctions of nanoribbons are formed. 

On a more fundamental level, the continuing great interest in the effect of reduced 
dimensionality, such as electron-electron interactions in reduced dimensions, provides 
motivation for studying quasi-one-dimensional systems. Graphene's unusual dispersion 
(massless Dirac fermions) and reduced density of states at the Fermi energy, for instance, 
suggest potential for novel effects. Of course, one should first understand the non-interacting 
system before turning to interactions. 

Graphene nanoribbons are closely analogous to electron waveguides patterned out of 
two dimensional electron gas (2DEG), usually in GaAs or other semiconductor systems 
ll32[33[Ml[35l[36l[37l[Ml[39l|40l|4Tl|^ However, there is an important difference in 
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how the confinement is achieved. While in 2DEG waveguides the electrons are trapped in 
the transverse direction of the waveguide by applying a potential by means of local gate 
electrodes, graphene nanoribbons are directly cut out of a larger graphene flake. This gives 
rise to different types of boundaries, depending on the direction in which the nanoribbons are 
cut with respect to the graphene lattice. If the longitudinal direction of the nanoribbon is along 
the direction of nearest neighbor carbon bonds, the resulting boundary is of "armchair" type, 
while cutting at 30° with respect to the nearest neighbor carbon bonds results in a "zigzag" 
boundary (see figure [TJ. It has been shown that the low energy properties of nanoribbons 
with boundaries other than these two are equivalent to those of zigzag nanoribbons |43|. On 
the experimental side, there has been recent progress in controlling the edges of graphene 
samples 14411451 . which is essential to enable physicists to probe the influence of edge details 
on transport properties. 

This paper is organized as follows: First we study one of the most simple systems beyond 
a straight nanoribbon with constant width, namely wide-narrow junctions, by which we mean 
two semi-infinite nanoribbons attached together to form a step. We calculate the conductance 
of such ribbons by numerically solving the tight binding model, and also obtain analytical 
results for the case of armchair boundaries. In the second part we investigate numerically the 
conductance of curved wires cut out of graphene. In this case the width of the nanoribbon is 
approximately constant, but the longitudinal direction with respect to the underlying graphene 
lattice and hence the transverse boundary conditions change locally. In contrast to systems 
with sharp kinks and abrupt changes in the direction, which have been investigated in earlier 
work |l7l[l3][l4l[T9l|20l|21 46|, we focus here on smooth bends. 

In both cases we find remarkable deviations from the conductance of 2DEG waveguides 
that are clear signatures of the sublattice and valley degrees of freedom in the effective 2D 
Dirac Hamiltionian describing graphene's low energy excitations, 

H^vJ "^p^ + ^yPy ^ )■ (1) 

Here the matrix structure is in valley space, (Tx/y are Pauli matrices in pseudo- or sublattice- 
spin space, pr^/y are the momentum operators, and vp ~ lO^m/s is the Fermi velocity. 
Alternatively, from a strictly lattice point of view, the deviations that we see are caused by 
the basis inherent in graphene's hexagonal lattice. 

For our numerical work, we use a nearest-neighbor tight binding model taking into 
account the 2p2-orbitals of the carbon atoms (|4] |42l and solve the transport problem using 
an adaptive recursive Green function method f48l to obtain the conductance G. Throughout 
the paper, lengths are given in units of the graphene lattice constant a which is V3 times the 
nearest-neighbor carbon-carbon length, while energies are in units of the nearest-neighbor 
hopping constant t = 2hvF/ (VSa) ~ 3 eV. 



2. Wide-narrow junctions: Changing the width of a nanoribbon 

The simplest way to form an interface within a nanoribbon is to change its width. In this 
section we investigate the conductance of infinite nanoribbons in which the width changes 
from wide to narrow, which then can be viewed as a junction between a wide semi-infinite 
nanoribbon and a narrow one. Figure[T]shows three examples of such junctions with armchair 
(ac) and zigzag (zz) type edges. We denote the width of the narrower wire by Wi and the 
width of the wider wire by W2 ■ A naive expectation for the dependence of G on the Fermi 
energy Ep is the step function G{Ef) ~ Ni {Ep)2e'^ /h where A^i is the number of occupied 
transverse channels in the narrow wire. This would be correct if there were no reflection at 




Figure 1. Wide-narrow junctions for different types of nanoribbons formed from a hexagonal 
lattice. The width of the nan'ower part is Wi while that of the wider part is W2 . The gray 
shaded sites denote infinite extension, (a) Abrupt junction between armchair nanoribbons. 
(b) Abrupt junction between zigzag nanoribbons. (c) Gradual junction between zigzag 
nanoribbons. 



the wide-narrow interface. Realistically, however, there is scattering from this interface, and 
so the steps in the conductance are not perfectly sharp. 

For usual 2DEGs modeled by either a square lattice of tight-binding sites or a continuum 
Schrodinger equation with quadratic dispersion, the detailed shape of G{Ep) has been studied 
previously. Szafer and Stone |36| calculated G{Ep) by matching the transverse modes of the 
two semi-infinite wires. The inset in figure|2](b) compares tight-binding results (using a square 
grid) with mode-matching results in this case for W2 — 2Wi in the one-mode regime of the 
narrow part. The agreement between the two is excellent. Note that the resulting conductance 
step is very steep. 



2.7. Armchair nanoribbons 

For armchair nanoribbons, the analysis proceeds in much the same way as for the usual 2DEG, 
square lattice case. At a fixed Fermi energy in the effective Dirac equation, the transverse 
wavefunctions for the various subbands are mutually orthogonal, as explained further in 



Appendix A Performing a matching procedure similar to that used in Ref. |36|, one calculates 



the conductance from the overlap of transverse wavefunctions on the two sides of the wide- 



narrow junction. A detailed derivation is presented in Appendix B 

Figure [2] shows the conductance resulting from the numerical solution of the matching 
equations at energies for which there is one propagating mode in the narrow part. In addition, 
the conductance obtained from tight-binding calculations for wide-narrow junctions between 
armchair nanoribbons is shown (using the hexagonal graphene lattice). Figure|2]shows G{Ef ) 
for different combinations of metallic and semiconducting nanoribbons (cf. [Appendix A) . 
The agreement between the two methods is extremely good: even the singularity associated 
with the subband threshold in the wider ribbon is reproduced in detail by the mode matching 
method, showing that the effective Dirac equation describes the system very well. 

In figure 121 we see immediately that G{Ep) for the armchair nanoribbon case differs 
greatly from the normal 2DEG G{Ep) [inset of figure |2](b)] : the rise from zero to unit 



Interfaces Within Gmphene Nanoribbons 

a) 





1 




0.8 - 








0.6 - 






CN 






4 - 


O 






0.2 - 




- 




1 




0.8 - 


J3 






0.6 - 


'a 




IN 






0.4 - 


o 






0.2 - 




- 




0.01 0.012 0.014 0.016 
C) 



0.005 0.01 0.015 0.02 0.025 
d) 




0.005 0.01 0.015 0.02 0.025 



0.012 



0.016 



0.02 



Ep [t] 



Figure 2. Conductance of wide-narrow junctions in armchair nanoribbons as a function 
of Fermi energy Ep : results from tight binding (solid lines) and mode matching (squares, 
obtained by solving equation )B.12) numerically). The energy window corresponds to the 
full one-mode regime of the narrow part. The behavior depends on whether the widths 
Wi^2 correspond to semiconducting or metallic armchair nanoribbons. (a) Semiconducting- 
semiconducting (Wi =99, W2 = 199). (b) Metallic-semiconducting (98, 199). (c) Metallic- 
metallic (98, 197). (d) Semiconducting-semiconducting (79, 109). Inset in (b): Conductance 
of a wide-narrow junction in a usual 2D electron gas: tight binding calculation (solid line) on 
a square lattice (Wi = 200 as, W2 = 400 as) and solution of matching procedure (circles, 
equation (2) of Ref. 1361 ). ts is the nearest neighbor hopping energy on the square lattice and 
as is its lattice spacing. 



conductance is much slower in graphene, taking at least half of the energy window and in some 
cases [see e.g. figure |2]( a)] not reaching the saturation value at all. For completely metallic 
nanoribbons, the lineshape is very different [panel (c)] and the conductance is suppressed at 
low Fermi energies (see also reference f49\). 

2.2. Zigzag nanoribbons 

For zigzag nanoribbons, the analysis does not proceed as simply as in the usual 2DEG 
or armchair nanoribbon cases: the transverse wavefunctions depend on the longitudinal 
momentum - similar to 2DEG wires with a magnetic field - and are not orthogonal at fixed 



Fermi energy (cf. Appendix A i. Because this orthogonality is used in the matching method of 
Appendix B[ we cannot apply it to the zigzag case. 

Figure [3](a) shows numerical tight binding results for G{Ef) in two different systems 
with zigzag edges: one with an abrupt change in width (red curve) and one with a gradual 
connection (blue curve), as depicted in figures [T](b) and (c), respectively. Note first that the 
conductance is close to its maximum value only in small windows of energy, as in the armchair 
nanoribbon case and in marked contrast to the usual 2DEG, square lattice situation. 

In the abrupt case, one sees pronounced antiresonances at the threshold energies for 
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Figure 3. (a) Conductance of wide-narrow junctions in zigzag nanoribbons as a function of 
Fermi energy (Wi ~ 19, W2 ~ 76). Curves for a structure with abrupt change in width 
[red line, depicted in figure 

Qb)] and one with a gradually changing width [blue dashed line, 
depicted in figure[T](c)] are compared to the number of propagating modes in the narrow wire 
(black line, i.e. the maximum possible conductance), (b) Probability density (color-coded in 
arbitrary units) of an electron entering the systein from the nan'ow region at Ep = 0.03 4. 
Only the density on the B sublattice is shown. A and B denote the sublattice type at the edges. 
The density decreases by a factor of about 20 from the B edge (red) to the A edge (blue). 



transverse modes in the wide nanoribbon. In order to see that this is due to the boundary 
conditions satisfied by the transverse modes in a zigzag nanoribbon, consider the following 
argument. As seen in Figures [T|(b) and (c), there is only one sublattice at each zigzag edge. 
In the effective Dirac equation one has a spinor with entries corresponding to the sublattices, 
thus the boundary condition is that one of the entries has to vanish at the edge while the other 
component is determined by the Dirac equation and is in general not zero at the boundary 
f9\. One finds from equation ( IA.24l i (e.g. from a graphical solution) that the higher Ep is 
above the threshold of a mode, the closer the transverse wavenumber gets to a multiple of 
tt/W and the closer the value of the spinor entry in question goes to zero. For our situation, 
then, the matching of a transverse mode in the narrow nanoribbon (which is already far 
above the threshold of the mode) with one in the wider nanoribbon is particularly bad at 
the threshold of the latter and gets better with increasing Fermi energy. This explains the 
observed antiresonances in G{Ef)- 

For the gradually widened junction, we insert another zigzag edge to interpolate between 
the wide and narrow nanoribbon [see figure[T](c)]. In this case, the modes of the two infinitely 
extended parts are not directly matched and thus the sharp antiresonances are not present. 
Note, however, the complete suppression of G at very low energies. In this regime there is 
only a single mode propagating in the wide nanoribbon as well as in the narrow one. This 
state is located mainly on the B sublattice close to the lower edge and on the A sublattice on 
the upper edge. Since the sublattice at the lower edge changes from A to B at the junction [cf. 
figure [T](c)], this state cannot be transmitted and the conductance is zero. This is confirmed 
by the intensity distribution plotted in figure|3](c). In the more realistic next-nearest-neighbor 
hopping model, the situation is the same for most of the single mode regime but changes 
for very low energies, when the so-called edge states are propagating In that regime, 

the two states are exponentially localized at the upper and lower edge, respectively, and are 
independent of each other. Thus, the one localized at the A edge transmits whereas the one 
localized at the B edge is blocked ll22l . 
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Figure 4. Schematics of curved graphene nanoribbons. zz denotes a nanoribbon with a zigzag 
edge as in figure[T](b), ac denotes an armchair edge, and "zz" denotes a zigzag-like boundary 
as explained in the text, (a) Parameters defining a "sidestep nanoribbon"; it is point symmetric 
about its center, (b) A structure with a single zigzag-armchair interface, (c) For 7 > 30°, 
there will be small regions with armchair edges (shaded red); these have a width Wac and 
behave as in |Appendix A| (d) For 7 < 30°, no armchair regions form; the curved nanoribbon 
is zigzag-like throughout. 

Summarizing the results for the wide-narrow junctions, we see that the behavior of 
graphene nanoribbons differs substantially from that of the familiar 2DEG situation. The 
matching at the graphene junctions is much less good, leading to a suppression of the 
conductance from the expected nearly step-Uke structure. 

3. Curved graphene nanoribbons 

Curved nanoribbons are defined by cutting smooth shapes out of an infinite graphene sheet. 
Since the graphene lattice is discrete, the resulting boundary is not perfectly smooth but will 
have edges of zigzag and armchair type in certain directions as well as some intermediate edge 
types. However, according to Akhmerov and Beenakker [43l, the intermediate boundary types 
behave basically like zigzag boundaries for low energies, and we thus call these boundaries 
"zigzag-like". 

In figurelHwe show schematically several of the curved nanoribbons studied. A "sidestep 
nanoribbon" consists of an infinitely extended horizontal zigzag ribbon of width W, followed 
by a curved piece with outer radius of curvature R2 and inner radius R1 — R2 — W, a second 
straight piece making an angle 7 with respect to the first one, a curve in the opposite direction, 
and finally followed by another infinitely extended zigzag nanoribbon. The details of the 
system's edge depend on 7: (1) If 7 = 30°, the middle straight piece has armchair edges. 
(2) If 7 > 30° the middle straight piece is zigzag-like with the dominating sublattice at the 
edges reversed from that for the two horizontal nanoribbons. In the curved part, there is a 
small region where the edges are locally of armchak type. If we denote the angle of the local 
longitudinal direction from the horizontal by 9, this happens at 6* = 30° [see figure|4](c)]. The 
inset in figure|5]shows the lattice structure of such a curved region. (3) Finally, if 7 < 30°, the 
middle straight piece also has zigzag-like edges, but the dominating sublattice at the edges is 
the same as for the horizontal ribbons. In this case, no local armchair region forms as 9 is 
always smaller than 30° [see figure|4](d)]. 

In these various cases, then, different interior interfaces are formed between zigzag and 
armchair nanoribbons. We will see that the type of interface is critical in determining the 
properties of the curved nanoribbons. In addition, the nature of the armchair nanoribbon - 
whether it is semiconducting or metallic - has a large effect on the conductance. Thus the 
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Figure 5. Conductance of sidestep nanoribbons as a function of Fermi energy with 7 = 60° 
and i?2 = 259 for two widths of the armchair region, VFac = 68.5 (solid red, corresponds to 
a semiconducting ribbon) and Wac = 68 (dashed blue, corresponds to a metallic ribbon). The 
solid black line shows the number of propagating transverse modes in the zigzag leads. Note 
that the internal interfaces between the zigzag and semiconducting armchair regions are much 
more reflective than for the metallic armchair case. Inset: The lattice structure of the first curve 
of a sidestep nanoribbon showing the armchair region formed at 6 = 30° . 



width of the armchair region W^c is an important parameter; according to equation ( IA.15I ) 
one has a metalUc nanoribbon if 4(1 + T4^ac/a)/3 G N and a semiconducting nanoribbon 
otherwise. 

Figure|5]shows the conductance of sidestep nanoribbons with 7 = 60°, for which a small 
armchair region is formed in each of the curved parts. When the width of this armchair region 
corresponds to a metallic nanoribbon, the conductance is basically 2e^//i - the maximum 
possible value - throughout the one-mode regime of the zigzag leads (W^c = 68, dashed blue 
line). In striking contrast, when the width is just a/2 larger (red line) the conductance is 
strongly suppressed. Resonance peaks result from Fabry-Perot behavior caused by scattering 
from the two armchair regions which define a "box" for the middle straight region. We find 
this behavior consistently for all sidestep wires in which armchair regions form that have a 
width corresponding to a semiconducting nanoribbon. 

Figure|6]shows the dependence on the angle 7 by plotting the conductance (G) averaged 
over all energies for which there is one propagating mode in the zigzag leads. For 7 < 30° 
there are no armchair regions in the curved parts of the structure, and the average conductance 
is very close to the maximum value in all cases studied. As soon as the critical angle of 30° 
is surpassed and small armchair pieces form in the curves, the conductance depends strongly 
on the exact value of W^c- If W^c corresponds to a metallic ribbon, (G) remains high and is 
rather independent of 7. On the other hand if W^c corresponds to a semiconducting ribbon, 
(G) suddenly drops by more than 80 percent and then remains approximately constant upon 
further increase of 7. The constancy of (G) in the respective regimes supports the statement 
of Ref. 1 43 1 that straight boundaries that are neither exactly of armchair nor exactly of zigzag 
type behave like zigzag boundaries. To summarize, if a curve in a zigzag nanoribbon causes 
two semiconducting armchair regions to appear, then a very effective barrier is formed which 
causes very high reflectivity. 
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Figure 6. Average conductance for two sidestep nanoribbons as a function of angle 7 
(R2 = 259). The average is taken over all Fermi energies in the one-mode regime of the 
zigzag leads. In one structure, the armchair region is metalHc (blue triangles, Wac = 68) 
while in the other it is semiconducting (red circles, Wa,c = 68.5). Note the sharp decrease in 
conductance in the semiconducting case when the armchair edges first form at 7 = 30° . The 
dotted lines are guides to the eye. 




Figure 7. Tight binding band structures of infinitely extended graphene nanoribbons. 
(a) Armchair nanoribbon with a width of 68.5 (same width as the local armchair piece 
forming in the structure of figure \5\ red curve), (b) Zigzag nanoribbon with a width of 
(39+ ^)\/3 ~ 68.4. The semiconducting energy gap in the armchair nanoribbon does not 
correspond to the energy region in figure|5]in which the conductance is suppressed. 



The simplest idea to explain this effect would be that at low energies there is by definition 
a gap in a semiconducting ribbon and since this means there are no propagating states in 
the local armchair region, one expects the conductance to be suppressed because electrons 
have to tunnel trough this region in order to be transmitted. However, this does not explain 
our findings: the energy range over which the conductance suppression occurs is much 
larger than the energy gap of the semiconducting region. In fact it is given by the energy 
range of the one-mode regime in the surrounding zigzag parts. To make this clear, we 
show the bandstructures of both a semiconducting armchair ribbon and a zigzag ribbon of 
approximately the same width in figure [T] (both nanoribbons are infinitely extended). One can 
clearly see that within the one-mode regime of the zigzag nanoribbon, in which the states are 
completely valley polarized, there can be several propagating modes in the semiconducting 
armchair nanoribbon, so the suppression of G must be of a different origin. 

Furthermore, it is not the bare zigzag-armchair junction that leads to suppressed 
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Figure 8. Single zigzag to amichair interface conductance of a smooth bend through 30°, 
as depicted in figure|4](b) (i?2 = 259). Both semiconducting (solid red line, VFac = 68.5) 
and metallic (dashed blue line, Wac = 68) armchair nanoribbons lead to good conductance. 
The solid black line shows the number of propagating transverse modes in the zigzag lead, 
coiTesponding to the maximum possible conductance (in the armchair lead at energies above 
the semiconducting gap there are always more or equally many inodes propagating). 




conductance, but rather it is necessary to have two zigzag pieces differing by an angle of more 
than 30° and being separated by a small armchair region. This can be seen in two stages. First, 
figure [8] shows the conductance of an infinitely extended zigzag nanoribbon connected to an 
infinitely extended armchair nanoribbon via a 30° curve, the structure shown schematically in 
figure|4](b). In the one-mode regime of the zigzag ribbon, the conductance is maximal for the 
case of a metalUc armchair ribbon. For a semiconducting armchair ribbon, the conductance 
is, of course, zero for energies below the gap, but it increases rapidly up to 2e^ /h for larger 
values of Ep. Thus, a single zigzag to semiconducting-armchair interface conducts well. 

For the second stage of the argument, consider a bend through 90° from an infinite 
zigzag lead to an armchair one, as depicted in figure|9] In contrast to the 30° zigzag-armchair 
connection just discussed, this one has three interfaces between zigzag and armchair regions. 
Figure [TO] shows the conductance of several 90° curved nanoribbons. As for the sidestep 
ribbons, the conductance is suppressed when a semiconducting armchair region is present in 
the curve. Note that the suppression is not due to the infinitely extended armchair lead, for 
which we chose a semiconducting nanoribbon in[T0](a) and a metallic one in[T0l(b). 

If one makes the natural assumption that the armchair region acts as a blocking barrier. 
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Figure 10. Conductance of 90° curved nanoribbons with either a semiconducting (solid red) 
or metallic (dashed blue) armchair region at 8 = 30° . (a) Semiconducting armchair lead: 
VFac = 69 (solid red) and Wac =69.5 (dashed blue), (b) Metallic armchair lead: Wac =88.5 
(solid red) and Wac = 89 (dashed blue). Black line: number of propagating modes in the 
zigzag lead. (iJ2 = 259.) 
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Figure 11. Average conductance as a function of the length added to the armchair region 
of a 90° curved nanoribbon. The structures are as in figure llOl (a): nanoribbons with a 
semiconducting (red) or metalHc (blue) armchair region at 8 = 30°. The average is over 
all Ep above the semiconducting gap of the armchair lead and in the one mode regime of the 
zigzag lead. The dotted lines are guides to the eye. 



one would expect the blocking to become more effective as the armchair region is lengthened. 
However, this is clearly not the case in the data shown in figure [TT] The system is a 90° 
curved nanoribbon in which the armchair region at 9 = 30° is lengthened by AL; we plot 
(G), the conductance averaged over all energies in the one-mode regime of the zigzag lead, 
as a function of AL. For a metallic armchair region in the curve, the conductance is roughly 
independent of AL, as expected. Surprisingly, for the semiconducting case, the conductance 
increases as a function of Ai.This establishes, then, that conductance suppression occurs 
when two zigzag-armchair interfaces occur in close spatial proximity. 

Our numerical results suggest that the evanescent modes in the armchair regions play an 
essential role. They are necessary, of course, in order to match the propagating zigzag mode 
to a solution in the armchair region. For short armchair pieces the evanescent modes from the 




Figure 12. Conductance of rings with armchair leads {R2 = 259). Red sohd: Ring with 
semiconducting armchair regions in both arms (Wc = 69 in the right arm and Wc = 68.5 
in the left). Blue dashed: Ring with metallic armchair regions in the light arm (Wc = 69.5) 
and semiconducting in the left (Wc = 69). Inset: Schematic of the ring structure; red shading 
indicates regions with armchair edges, as in figure|4] 



two interfaces overlap. We conjecture that these evanescent modes are mutually incompatible 
in the semiconducting case, destroying the possibility of matching on both sides at the same 
time, while they are compatible for metallic armchair regions. If one has a long armchair 
piece, the evanescent modes decay leading to independent matching at the two ends. 

4. Conclusions 

We have shown in a variety of examples that interfaces within graphene nanoribbons 
can strongly affect their conductance, much more so than in the familiar 2DEG electron 
waveguides and wires. First, for wide-narrow junctions, our main results are figures |2] and 
[3] For both armchair and zigzag nanoribbons, changes in width act as a substantial source of 
scattering, reducing the conductance. Second, for curved nanoribbons, our main results are 
figures |6] [8] andfTTl There is a strong reduction in conductance when a curve joining two 
zigzag regions contains a semiconducting armchair region. 

The effect of such internal interfaces will certainly be felt in more complex structures 
as well. As an example, consider rings for studying the modulation of the conductance 
by a magnetic field through the Aharonov-Bohm effect |50|. Figure [T2l shows such a ring 
schematically together with its conductance in two cases. As for the curved nanoribbons, 
when semiconducting armchair regions occur in the curved part of the structure, the 
conductance is substantially reduced. 

In considering experimental manifestations of internal interfaces, disorder, and in 
particular the edge disorder which has received attention recently 1441 l45l ISTl l52l l53l 1541 . 
may be important. The effects we observe in our calculations will most likely also be present 
in structures with disordered edges, provided the disorder is not too strong. Consider, for 
example, the suppression of the conductance in curved wires. In the inset of figure|5]as well 
as in figure |9l one can see that between the armchair and the zigzag regions, the edges are 
not perfect. We believe that when the edge disorder is weak enough to allow for pieces with 
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armchair edges, the suppression should still be present. 

The underlying reason for the impact of internal interfaces can be viewed in two ways. 
From the lattice point of view, it arises from the additional complexity of the hexagonal lattice 
with its basis compared to the standard square lattice. Equivalently, from the continuum 
point of view, it arises from the extra degrees of freedom inherent in the Dirac-like equation 
governing graphene - those of the sublattice and valley pseudospins. As development of 
graphene nanostructures accelerates, the impact of internal interfaces should be taken into 
account when considering future carbon nanoelectronic schemes. 
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Appendix A. Wavefunctions of graphene nanoribbons in the Dirac equation 



We calculate the eigenfunctions of graphene nanoribbons withing the effective Dirac model. 
This has been done by Brey and Fertig in |9| and Peres, et al in ifTOll . The effective Dirac 
equation taking into account contributions from both valleys is given by ID 



(A.l) 



x = 



a) 




x = W 




y = 



: Sublattice A 
: Sublattice B 



Figure Al. Infinitely extended graphene nanoribbons. (a) Armchair ribbon along the y- 
direction. The outermost rows of atoms are at x = a/2 and W — a/2 respectively. Hence, 
the width of the ribbon is given by W = W — a. The boundary condition however is, that 
the wavefunction is zero at a; = and x = W respectively, (b) Zigzag ribbon along the x- 
direction. Here the width of the ribbon isW = W — 2a/ V3. Since first row of missing atoms 
at each side is only on one sublattice, the boundary conditions requires only the corresponding 
part of the wavefunction to vanish. 
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with 

\^ -Crj:Pa: + GyPy J 

and 

$(r) = [<^K{r),<^K'{r)f = [^^(r), $s(r-), ^^iM, $'s(r)]^ . (A.3) 

Here $x and are spinors with two components corresponding to contributions from the 
two different valleys K and K' respectively, ^a/b and ^'^^g are scalar wavefunctions, where 
the subscripts A and B stand for the two sublattices (see figure IATT i. The total wavefunction 
containing the fast oscillations from the iiT-points is then 



Appendix A.l. Armchair nanoribbons 

We consider an armchair nanoribbon which is infinitely extended along the y-direction [see 



figure |AT]( a)]. Using the Bloch ansatz 

$(r) = e:''yy(t){x) (A.5) 
and the Dirac equation (lA.U . one obtains 

- i{ky + dx)(l}B{x) ^ e(l)A{x) (A.6) 

i(ky - dx)4'A{x) = e(l>B{x) (A.7) 

-t{ky - d^)^'B{x) ^ e^'Aix) (A.8) 

t{ky + d^)cj,'^{x) ^ e^'six) (A.9) 
and, by applying the Hamiltonian twice, 

ikl~dl)(t>ix)^e^^{x) (A. 10) 



with e = E/{hvF)- According to figure |AT]( a), the correct boundary condition ||9j for an 
armchair nanoribbon is V'(^) ~ for x = and x ~ W. (For the connection between 
the nanoribbon width W used previously and W, see the caption of figure IaTI ) The ansatz 
(jjBix) = Ae*«"^ + Se-*«"^, <j)'s{x) = Ce"?"^ + De""?"^ solves both the B sublattice parts 
of equation (lA.lOl i with = fc^ + qf^ and the boundary condition, if we require 

ni: „ 
(?„ = -: K with n e Z (A. 1 1) 

where K = 47r/(3a). We find that B = C = and A = —D. Using equations iA.6\ and 
(IA.8b to determine <i>yi(a;) and from $5(2;) and $^(a;), we thus find that, up to a 

normalization factor, the wavefunctions are 

</.(x) ~ [(g„ - iky)e''^-ye, e''"^ -(g„ - ifc,)e-''"7e, -e-^^-^]^ , (A.12) 

^(r.)^e''=«^sin[(g„ + if)x][(g„-ifcy)/e, 1]^. (A.13) 

The wavefunction il'{r) is, up to the spinor part, very similar to that of a 2DEG 
waveguide: the width of the ribbon is a multiple of half the transverse wavelength. However, 
here the transverse wavelength is of order the lattice constant, not the system's width, since n 
is of order W/a for the energetically lowest lying modes. Nevertheless, the wavefunctions for 
different transverse quantum numbers n are orthogonal at a fixed Fermi energy. Note that for 
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evanescent modes we just have to consider imaginary wavenumbers ky = iKy and equations 
jATTT i, ( |Al2] i, and jXlll i still hold. 
The energy for this solution is 



E = ±hvF^kl + ql. (A.14) 

Therefore one has a metallic spectrum if there is a state with (/„ = 0. From equation (lA.lll i it 
follows immediately that this is the case whenever 

4T4^ 

-— eN. (A.15) 

6 a 

Appendix A.2. Zigzag nanoribbons 

For a zigzag nanoribbon along the x-direction [see figure|AT|(b)], the Bloch ansatz is 

$(r) = e''''="=0(?;) , (A. 16) 

the Dirac equation becomes 

(k, - dy)(t>Biy) ^ e<f>Aiy) (A.17) 

[k, + dy)<f>A{y) = eMy) (A.18) 

-{k, + dy)<j>'siy) ^ e<j>'Ay) (A19) 

-{k, - dy)^'Ay) ^ ^^'siy) , (A.20) 

and one has 

[kl - dl)m = e^ct^iy) . (A21) 

The boundary condition for a zigzag ribbon differs from that for an armchair ribbon in 
that the wavefunction has to vanish on only one sublattice at each edge [§]: ipAix, y = 0) = 
ipsix, y — W) = 0. With the following ansatz for ^a{x) and $^(a;), 

(y) = Ae'^y + Be-^'-y (jy'^ (y) = Ce^^ + De-''^ , (A.22) 

( IA.2II 1 yields = k'^ + z'^, and the boundary condition requires A — —B and C ~ —D. Thus 
the valleys completely decouple for zigzag nanoribbons, and equations ( IA.I8I 1 and ( IA.20l i 
yield 

"^K/K- ^ [sinlzy), {r fc^r sin(2;y) ^ z cos.{zy)} / ef (A.23) 

where r = + 1 for the K and t = — 1 for the K' valley. The boundary condition for the B 
parts of the wavefunction provides an equation that determines the allowed values for z, 

k^ = -Tz/i{m{zW) . (A.24) 

Thus the transverse quantum number is coupled to the longitudinal momentum, as in 2DEG 
waveguides in the presence of a magnetic field. In order to write equation ( IA.23b in a 
symmetric way, we square the quantization condition ( |A.24| l and use the relation fc^ = — 
to obtain 

= z'^/siT?{zW). (A.25) 
Using (IA.24I I and ( |A.25l l in equation ( |A.23l l leads to 

^K/K' - [sin(z2/), s{z, e) sin{z(W^ - y)}f (A.26) 

with s(z, e) — sign[e2;/ sin(zM^)]. ^From this symmetric expression, one clearly sees that the 
total weight on each sublattice is the same. 



15 




The transcendental equation ( IA.24l i has real solutions z e M only for |e| > 1/W. These 
states correspond to bulk states: they are extended over the whole width of the ribbon. For 
|e| < 1/W there are only imaginary solutions iz € M, corresponding to the so-called edge 
states l6]|9l, which are exponentially localized at the edges and live predominantly on one 
sublattice at each side, as can be seen from equation ( IA.23I ). For the special case 2 = 
corresponding to |e| = 1/W equation ( IA.26I I results in 

lim ^K/K' - [y, -sign(e)(y - W)f , (A.27) 

i. e. a linear profile of the transverse wavefunction. Figure (IA2I ) shows the profile of several 
transverse zigzag modes. 



Appendix B. Mode matching for wide-narrow junctions with armchair edges 

We derive a set of analytic equations that determine the transmission amplitudes for wide- 
narrow junctions with armchair edges as introduced in section [TT] We label the transverse 
modes in the narrow part of the system by f^{x) and those in the wide part by x^i^)- The 
± stands for propagation in positive and negative y-direction respectively. Furthermore we 
use latin subscripts n and m for ip and greek subscripts i' and oj for x- Then we know from 
[Appendix A| that 

iftix) ^^sinimrx/Wi) [(q„Tifc^)A, l]^ (B.l) 
V Wi 

xtix) = ^sin{yTix/W2) [{q,Tikl)/t, l]^ . (B.2) 

Here, we define the ky^" to lie always on the positive real axes for propagating states and on 
the positive imaginary axes for evanescent states 
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The full scattering wavefunction for an electron incident from the wide side in mode ui 

is 

y<0: Mx,y) = xi(a^)e'''«^ + 5I^-x;(2:)e-''=«^ (B.4) 



y 



where the sums run over all modes, both propagating and evanescent. Matching the two parts 
at the junction, defined to be y = 0, we obtain 

u n 

We can extract the scattering amplitudes by projecting this equation first on the wide side and 



and 



then on the narrow side. First, multiplying the B-part of this equation by 
integrating from to W2 yields 

-S,^ + ^ 2tnujb-+ (B.7) 



for which we used dx Xv' si^) I si^) — h^i^^' the definition 



W2 

htt ■■= J dx 





Since f^si^) vanishes for x > Wi, one can replace the upper limit of integration W2 by 
Wi. 

Second, we project equation iB.6i onto modes of the narrow lead. Multiplying by 
[ip^, (a;)] and integrating from to Wi yields 

where we have introduced the definitions (note the spinor inner product) 

d±±:=/ dx{^i{x)\ xt{^) (B.IO) 
Jo 

and have again used orthogonality of the transverse wavefunctions, now in the form 

dx{^l(x)\ ^+(:c) = — (|(7„ + zA;;f + (B.ll) 
Combining equations ( IB.7I I and ( IB.9I ), we obtain 



^ ( ^ 2d+-6;+ - ^ (|g„ + z/c^p + <5„„, ] t„,^ = rf+J - d++ (B.12) 
which can be written as a matrix equation in the form 

M„m traui = C„ . (B.13) 

m 

This can be solved for the by introducing large enough cut-offs for to and v and then 
inverting the now finite matrix M. 
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The total transmission for a particle incident in mode to from the wide side is given by 



n prop. 



n prop. 



7 



\tr. 



E 

n prop. 



Ky 



\tr. 



(B.14) 



Finally, the conductance of the system is cormected to the transmission via Landauer's formula 

2e2 



OJ prop. 

In these last two equations the sums run over propagating modes only 



(B.15) 
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